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ABSTRACT 

The  linear  stability  of  compressible  plane  Couette  flow  is  investigated.  The  correct  and 
proper  basic  velocity  and  temperature  distributions  are  perturbed  by  a  small  amplitude  nor¬ 
mal  mode  disturbance.  The  full  small  amplitude  disturbance  equations  are  solved  numeri¬ 
cally  at  finite  Reynolds  numbers,  and  the  inviscid  limit  of  these  equations  is  then  investigated 
in  seme  detail.  It  is  found  that  instability  can  occur,  although  the  stability  characteristics  of 
the  flow  are  quite  different  from  unbounded  flows.  The  effects  of  viscosity  are  also  calculated, 
asymptotically,  and  shown  to  have  a  stabilizing  role  in  all  the  cases  investigated.  Exceptional 
regimes  to  the  problem  occur  when  the  wavespeed  of  the  disturbances  approaches  the  ve¬ 
locity  of  either  of  the  walls,  and  these  regimes  are  also  analyzed  in  some  detail.  Finally,  the 
effect  of  imposing  radiation-type  boundary  conditions  on  the  upper  (moving)  wall  (in  place 
of  impermeability)  is  investigated,  and  shown  to  yield  results  common  to  both  bounded  and 
unbounded  flows. 
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1.  Introduction 


Incompressible  plane  Couette  flow  possesses,  perhaps,  the  simplest  exact  solution  of  the 
Navier-Stokes  equations,  and,  (probably  as  a  consequence)  the  study  of  the  stability  of  this 
flow  has  been  the  subject  of  considerable  attention  over  the  years.  Numerical  studies  of  the 
linear  stability  problem  have  been  carried  out  by  Grohne  (1954),  Gallagher  and  Mercer  (1962, 
1964),  Deardorff  (1963),  Davey  (1973),  and  Gallagher  (1974);  however,  all  these  studies  found 
no  evidence  of  instability.  A  number  of  analytic  studies  have  also  been  carried  out  on  this 
problem.  Wasow  (1953)  showed  that  the  flow  is  stable  at  all  streamwise  numbers  (a)  if  the 
Reynolds  number  (Re)  is  sufficiently  large.  The  stability  of  the  flow  at  low  Reynolds  number 
was  demonstrated  by  Synge  (1938).  Dikii  (1964)  proved  that  all  modes  with  wavespeeds  equal 
to  the  average  of  the  wall  velocities  were  stable  and  indeed  that  the  imaginary  component 
of  the  wavespeed  Im{c}  <  — a/Re.  The  first  general  proof  of  stability  appears  to  be  due  to 
Romanov  (1973).  He  showed  that  all  normal  modes  of  the  linear  problem  are  damped  for 
a  >  0,  Re  >  0.  Exact  solutions  of  the  Orr-Sommerfeld  equations  have  been  obtained  by 
Reid  (1979). 

Another  issue  that  has  been  studied  is  the  question  of  a  continuous  spectrum.  Case 
(1960,  1961)  showed  that  the  time  dependent  inviscid  problem  has  a  continuous  spectrum 
which  decays  in  time  as  1/t,  this  spectrum  arising  as  a  direct  consequence  of  the  singularity 
of  the  inviscid  equations  at  the  critical  layer.  Shivamoggi  (1982)  presented  an  example  of  a 
continuous  spectrum  which  decayed  as  an  exponential  in  time. 

A  detailed  analysis  of  modes  for  the  large  Reynolds  number  limit  for  general  mean  flows  is 
given  by  Morawetz  (1952),  which  has  implications  for  plane  Couette  flow.  She  classified  the 
modes  into  three  sets.  In  the  first  set  are  eigenvalues  which  approach  the  inviscid  eigenvalues 
in  the  limit  of  infinite  Reynolds  number.  Since  there  are  no  discrete  inviscid  eigenvalues  in 
the  case  of  plane  Couette  flow,  this  first  set  is  empty.  In  the  second  set,  there  ;s  an  infinity 
of  eigensolutions  unrelated  to  the  inviscid  problem,  satisfying 

Ic-Cnl  <  A(aRe)"2,  (1.1) 


where  A  is  a  constant,  and  Cn  is  a  root  of 
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where  n  is  any  integer,  Uo(y)  is  the  velocity  profile,  a  is  the  wavenumber,  and  the  flow 
extends  from  y  =  y\  to  y  —  These  modes  are  always  stable  (if  Re  is  sufficiently  large). 
In  the  third  class  are  eigensolutions  of  the  viscous  problem  that,  as  aRe  — ►  oo,  approach  a 
finite  V-shaped  strip  in  c-space  defined  by  one  branch  of 

Rl  r[i(Uo(y)-c)}tdy  =  0  (1.3) 
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and  one  branch  of 

Rl  I  [i{U0(y)  -  c)]%dy  =  0, 
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where  Uo(yc)  =  Rl{c}.  These  modes  are  either  neutral  or  stable. 
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However,  there  does  exist  a  dichotomy  between  the  theoretical  and  computational  results 
described  above  and  experimental  results  at  large  Reynolds  numbers,  in  which  instabilities 
are  certainly  observed  (Taylor  193d,  Reichardt  1956,  Robertson  1959).  A  number  of  attempts 
have  been  made  to  explain  this  through  the  inclusion  of  nonlinear  terms.  Investigations  of 
this  type  include  the  work  of  Watson  (1960),  Eckhaus  (1965),  Hains  (1967),  Kuwahara 
(1967),  Reynolds  and  Potter  (1967),  Ellingsen,  Gjevik,  and  Palm  (1970),  and  Lessen  and 
Cheifetz  (1975).  These  studies  used  a  variety  of  techniques,  including  varying  degrees  of 
mathematical  rigor,  and  led  to  a  number  of  conclusions  (some  partly  contradictory),  although 
on  balance,  the  evidence  was  that  finite  amplitude  effects  could,  indeed,  destabilize  the  flow. 
Additionally,  Romanov  (1973)  considered  the  nonlinear  initial-value  problem  and  showed 
that  there  is  a  unique  solution  which  is  asymptotically  stable  if  the  norm  of  the  initial 
disturbances  in  Sobolev  space  is  sufficiently  small. 

The  linear  stability  of  compressible  flows  is  considerably  less  well  understood  than  corre¬ 
sponding  incompressible  flows.  Most  of  the  work  (which  has  for  the  most  part  been  based  on 
the  parallel  mean  flow  assumption)  has  been  with  regard  to  boundary-layer  flows  (e.g.,  Lees 
and  Lin  1946,  Mack  1963,  1965a, b,  1969,  1984,  1987,  1990)  and  more  recently  to  jets  and 
shear  layers  (e.g.,  Tam  and  Hu  1989a,b,  Greenough  et  al.  1988,  Papageorgiou  1990,  Jackson 
and  Grosch  1989,  Mack  1990). 

The  first  (of  a  number)  of  distinguishing  features  of  the  stability  of  compressible  flows 
was  found  by  Lees  n."d  Lin  (1946).  This  relates  to  the  replacement  of  the  (classical)  inflec¬ 
tion  point  condition  In  the  streamwise  velocity  profile  (for  the  existence  of  neutral  inviscid 
modes),  by  the  generalized  inflection  point  (G.I.P.)  condition  which  involves  the  mean  den¬ 
sity  distribution  also.  This  condition  relates  to  the  existence  of  neutral  “inviscid”  modes  (i.e., 
modes  of  wavelength  comp?.raf  i*  fo  the  characteristic  scales  of  the  mean  flow,  for  example 
the  displacement  thickness  in  the  case  of  bound  ary- layer  flows,  or  lateral  extent  in  the  case 
of  confined  flows). 

Another  effect  of  compressibility  is  the  frequent  existence  of  multiple  modes  of  instability 
in  supersonic  flows.  In  addition  to  the  modes  r«J  :+?d  to  the  G.I.P.,  a  further  class  of  modes 
may  exist;  these  are  linked  to  the  so-called  supersonic  neut.al  modes  (as  opposed  to  subsonic 
neutral  modes  associated  with  G.I.P. ’s).  Such  neutral  modes  are  readily  distinguishable, 
because  of  the  oscillatory  nature  of  the  eigensolutions  in  the  free  stream,  while  subsonic 
neutral  modes  exhibit  exponential  decay  of  their  eigensolutions  in  the  free  stream.  Mack 
(1963,  1965a, b,  1969,  1984,  1987,  1990)  has  shown  how  such  modes  are  extremely  important 
in  compressible  boundary-layer  flows,  particularly  in  the  case  of  boundary  layers  on  cooled 
walls,  where  these  modes  are  often  associated  with  the  largest  amplification  rates.  Such 
modes  have  also  been  shown  to  exist  and  be  important  in  other  classes  of  flow  by  Tam 
and  Hu  (1989),  Greenough  (1989),  Zhuang,  Kubota,  and  Dimotakis  (1990)  (confined  two- 
dimensional  supersonic  mixing  layers),  Macaraeg  (1990),  and  by  Jackson  and  Grosch  (1989) 
(compressible  mixing  layers).  However,  to  date,  a  simple  mechanism  that  explains  why  and 
precisely  under  what  conditions  these  types  of  mode  exist  has  yet  to  be  propounded.  Some 
general  remarks  concerning  stability  at  large  Reynolds  numbers  have  been  made,  however, 
by  Morawetz  (1959),  in  particular,  the  presence  of  inner  viscous  regions  and  their  nature 
was  vigorously  established  in  both  the  incompressible  and  compressible  cases.  In  the  case 
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of  stable  disturbances,  with  finite  damping,  it  was  shown  that  a  finite  viscous  region  in 
t'.?  interior  of  the  fluid,  no  matter  hov  large  the  Reynolds  number,  exists.  In  the  case 
of  unstable  disturbances,  with  finite  amplification,  there  is  no  inner  viscous  region  if  the 
Reynolds  number  of  the  mean  flow  is  large  enough.  In  the  case  of  neutral  disturbances  (for 
constant  Prandtl  number)  it  was  also  rigorously  shown  that  if  a  viscous  region  does  exist, 
its  width  shrinks  to  zero  with  an  increase  in  Reynolds  number. 

The  object  of  this  study  is  to  analyze  some  of  the  characteristics  of  compressible  plane 
Couette  flow.  Although  analytic  expressions  for  the  mean  streamwise  velocity  Uo(y)  and 
temperature  To(y)  are  not  available  for  general  viscosity  laws.  To(y)  may  be  expressed  as 
a  second-order  polynomial  in  Uo(y ),  and  hence  it  is  relatively  straightforward  to  generate 
mean  flow  piofiles,  under  different  conditions.  Some  work  has  previously  been  carried  out 
investigating  the  compressible  stability  of  plane  Couette  flow.  Glatzel  (1989)  has  considered 
this  problem,  and  has  carefully  studied  the  least  stable  eigenmodes.  He  also  derived  critical 
Reynolds  numbers  for  the  instabilities  that  were  found.  However,  this  work  is  for  the  special 
case  of  constant  viscosities,  density  and  pressure  (implying  constant  temperature  also),  which 
lead  to  a  fourth-order  system  (rather  than  a  sixth-order  system).  Girard  (1988)  considered 
the  same  problem,  and  also  assumed  constant  viscosities  (although  he  did  allow  for  variations 
in  density  and  temperature);  instaoilities  were  again  found  to  occur.  In  this  paper,  we  choose 
to  consider  the  stability  of  the  mean  flow  profile  described  above,  which  is  a  correct  and 
proper  solution  of  the  full  compressible  Navier-Stokes  equations.  We  show  that  the  details 
of  the  mean  flow  profile  have  a  profound  effect  on  the  stability  of  the  flow. 

The  layout  of  the  paper  is  as  follows.  In  Section  2  we  formulate  the  problem,  stating 
our  fundamental  assumptions  and  equations  of  motion.  In  Section  3  we  derive  the  equations 
of  motion  for  the  basic  flow,  in  which  the  problem  is  reduced  to  that  of  a  straightforward, 
although  nonlinear,  first  order  system  that  may  be  solved  by  means  of  standard  numerical 
means.  In  Section  4  we  derive  the  (full)  small  amplitude  disturbance  equations,  neglecting 
only  terms  in  perturbation  amplitude  squared,  and  we  describe  a  numerical  scheme  to  treat 
this  system,  together  with  a  number  of  numerical  results.  In  Section  5  we  consider  the 
inviscid  limit  of  these  equations;  we  show  how  the  so-called  “generalized  inflection  point”  is 
relevant  in  this  context,  and  under  what  circumstances  we  can  expect  such  a  point  to  occur 
with  our  basic  profile. 

Numerical  results  for  the  inviscid  problem  are  presented  in  Section  6,  and,  most  l  npor- 
tantly,  it  is  shown  that  unstable  modes  are  possible,  although  according  to  inviscid  theory, 
there  are  many  regions  where  many  modes  are  neutrally  stable.  These  neutral  modes  are 
investigated  further  in  Section  7,  in  which  the  effects  of  viscosity  are  included,  and  shown 
to  always  have  a  stabilizing  role  for  high  Reynolds  numbers.  However,  this  study  also  raises 
important  questions  regarding  the  applicability  of  our  results,  in  particular  in  regions  of 
changeover  from  the  neutral  to  non-neutral  state.  This  region  is  investigated  in  some  detail 
in  Section  8.  Because  of  the  apparent  discrepancy  in  results  between  bounded  and  un¬ 
bounded  stability  analyses,  in  Section  9  we  consider  a  change  in  boundary  conditions  on  the 
upper  (moving)  wall,  from  one  of  impermeability  to  one  of  radiation.  Finally,  in  Section  10 
we  present  some  conclusions. 
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2.  Formulation 


We  assume  that  we  have  a  compressible,  Newtonian,  perfect  flaid  between  two  infinite 
parallel  planes  defined  by  y*  =  0  and  y*  =  h.  The  x*  axis  is  taken  to  lie  in  the  plane  of  the 
lower  wall.  We  take  the  fluid  to  have  density  p',  viscosity  /if*,  second  coefficient  of  viscosity 
C  (which  may  be  taken  to  be  zero  for  a  mono  tonic  gas).  The  upper  wall  has  velocity  (/<*, 
(a  subscript  oo  refers  to  unperturbed  conditions  on  the  upper  wall),  taken  parallel  to  the 
plane  of  the  wall,  while  the  lower  wall  is  at  rest.  Although  the  basic  flow  will  be  taken  to 
be  one- dimensional,  later  we  shall  consider  two-dimensional  perturbations  of  this  flow.  The 
velocity  components  are  taken  to  be  u*  =  (it*,  vm)  in  the  x*  and  y*  directions  respectively,  and 
the  pressure  and  temperature  are  written  as  p*  and  T *  respectively.  We  non-dimensionalize 
velocities  with  respect  to  Uoo,  lengthscales  with  respect  to  h,  dens' ty  with  respect  to  p^, 
viscosities  with  respect  to  p^,  temperature  with  respect  to  and  pressure  with  respect 
to  Pto-^*^oo,  where  the  gas  constant  R*  =  Cp  —  Cv,  and  Cp  and  Cv  are  the  specific  heats 
at  constant  pressure  and  volume  respectively;  non-dimensional  quantities  are  denoted  using 
the  same  notation  as  corresponding  dimensional  quantities,  except  without  the  superscript 
asterisk. 

The  continuity  equation  may  then  be  written 


<21) 

Here,  and  throughout  the  paper,  we  assume  that  the  appropriate  dimensional  timescale  is 

Oih/Uco). 

The  momentum  equations  are  written 
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Here  the  Reynolds  number  Re  is  defined  by 
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the  Mach  number  M ^  by  _ 

=  UjJWri,,  (2.5) 

and  the  ratio  of  specific  heats  by  7. 
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The  non-dimensional  form  of  the  energy  equation  used  is 


dT  dT  3T]  ( 7  — 1\ 


dp  dp  dp  1  [  3 
dt  Udx  +  Vdy  Re|dy 
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Here  a  is  the  Prandtl  number 

<r  =  p-Cp/K', 

where  K*  is  the  coefficient  of  heat  conductivity.  The  equation  of  state  is  simply 


(2.6) 

(2.7) 


p  =  pT.  (2.8) 

We  assume  that  viscosity  depends  solely  on  temperature,  and  in  particular  we  assume  Suther¬ 
land’s  law 

/•  =  r^(i±£).  (2.9) 

where  C  is  a  constant.  Finally,  A  =  £  —  2/3  p  and  the  Stokes  assumption  £  =  0  is  assumed 
throughout  the  paper. 

In  the  following  section  we  consider  the  basic  flow  which  we  expect  to  depend  on  y  only. 


3.  Compressible  Couette  Flow 

We  seek  a  solution  of  (2.1)  -  (2.6)  which  is  dependent  on  y  only,  together  with  a  constant 
mean  pressure.  By  continuity  we  must  have  v  =  0.  We  then  seek  a  solution  of  the  form 


«  =  Uo{y),  T  =  T0{y),  p  =  p0(y). 

We  then  have 

{P0Ucy)y  =  0, 

subject  to 

t/o(0)  =  0,  U0{  1)  =  1, 

To( o)  =  rw>  r0(i)  =  i. 

It  follows  immediately  from  (3.2)  that  the  shear  stress  r  is  a  constant  through  the  profile, 
i.e., 

r  =  poUoy  =  constant.  (3.5) 

The  energy  equation  may  then  be  written  as 


(3.1) 

(3.2) 

(3.3) 

(3.4) 


—Ton  +  (7  -  mLUor 

C  Jy 


=  0. 


(3.6) 
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Defining  the  recovery  temperature  (i.e.,  the  wall  temperature  with  adiabatic  conditions) 


and  the  recovery  factor 


(3.6)  may  be  integrated  and  w  'itten  in  the  form 


(3.7) 

(3.8) 

(3.9) 


Equations  (3.5)  and  (3.9),  together  with  an  appropriate  viscosity  law  (we  used  (2.9))  and 
boundary  conditions  (3.4)  then  completely  determine  the  problem;  r  is  unknown  a  priori, 
and  must  therefore  be  determined  as  part  of  an  iterative  process.  The  mean  pressure  profile 
is  constant. 


4.  The  Small- Amplitude  Disturbance  Equations 

Here  we  take  the  solution  to  be  that  of  Couette  flow  (as  discussed  in  the  previous  section) 
together  with  a  small  amplitude  perturbation.  More  specifically,  we  write 

u  =  U0(y)  +  6u(y)E  +  O(62), 
v  =  6av(y)E  +  0(62), 

P  =  1  +  6p(y)E  +  0(6?), 

P  =  Po(y)  +  6p(y)E  +  0(<52),  (4.1) 

P  =  Po  (y)  +  6jx(y)E  +  0(6% 

C  =  Co  (y)  +  6C(y)E  +  0(6% 

T  =  T0(y)  +  6T(y)E  +  0(6% 

where 

E  =  exp  [m(x  —  d)] ,  (4.2) 

and  6  is  the  (small)  amplitude  of  the  perturbation.  If  we  take  terms  0(6)  in  (2.1)  -  (2.3), 
(2.6)  we  obtain  the  following  full,  small  disturbance  equations: 

-icp  -f  ip0u  +  iUop  +  vpoy  +  poVy  =  0,  (4.3) 


Po  \-iacu  +  iaUou  +  avUoy\  + 


j-2 p0a2u  +  A0  [-a2u  +  iauy]  +  p0y  [uv  +  iav ]  +  p0  [uyy  +  iavy]  +  pyU0y  +  pUoyy}, 

6  (4-4) 


Po  -ia2v  +  ia2Uov  + 


~  {2poyVy  +  2povyy  4  ^Oy  H  +  vy)  +  A0  (iauy  +  vyy)  +  po  (iotUy  -  a2#)  +  iotpUoy }, 

(4.5) 
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(4.6) 


p0  \j-iacT  +  iaUoT  +  auT0y]  —  ~  ^  [—iacp  +  iaUop]  = 

— \jlQ>Tvy  +  flOyTy  +  pToyy  +  fiyTQy  ~  0? floT  j 
+^— ^  [pol'oyUy  +  p(U0y)2  +  2fi0iaU0yV }  . 


Re 

The  perturbation  equation  of  state  is  then 


P  =  Top+  jT- 


(4.7) 


No-slip  and  impermeability  conditions  are  applied  to  the  velocity  perturbations  at  both  walls, 
namely 

5(0)  =  fi(l)  =  »(0)  =  5(1)  =  0.  (4.8) 

The  temperature  boundary  conditions  are  T(l)  =  0,  while  on  the  lower  wall, 


for  a  heated/cooled  surface  and 


dT 

dy 


T  (0)  =  0 


=  0  on  y  — 


(4.9) 


(4.10) 


for  an  insulated  wall. 

We  now  present  results  based  on  the  full  set  of  viscous  compressible  linearized  equations, 
i.e.  (4.3)-(4.7)  above. 

The  stability  results  are  obtained  from  a  spectral  temporal  linear  stability  code  written 
by  Herbert  (1990).  The  code  assumes  a  global  representation  of  all  variables  in  appropriate 
basis  functions  which  can  vary  from  variable  to  variable.  The  basis  functions  for  the  velocity 
components  satisfy  the  imposed  boundary  conditions.  Velocity  perturbations  are  set  to  zero 
at  y  =  0  and  y  =  1.  Thus,  the  velocity  basis  functions  are  linear  combinations  of  Chebyshev 
polynomials.  In  terms  of  the  Chebyshev  polynomials  Tn{y),  these  basis  functions  are  defined 
as 

Un(y)  =  Tn+2(y)  -  Tnk).  (4.11) 

Boundary  conditions  for  the  density  and  temperature  perturbations  waves  are  of  more  general 
nature,  so  these  variables  are  represented  as  a  standard  Chebyshev  series.  In  this  section, 
only  adiabatic  results  are  shown.  Thus,  the  y  derivative  of  the  temperature  perturbation  is 
zero  at  the  lower  wall.  The  upper  wall  is  insulated;  therefore  the  temperature  perturbation 
is  set  to  zero.  The  density  at  the  walls  is  obtained  by  integration  of  the  continuity  equation. 

The  numerical  algorithm  departs  from  standard  collocation  schemes  which  operate  in 
physical  space.  Taking  the  perturbation  density  ( p )  as  an  example,  let 


P(y)  =  ^2anTn{y)j 


(4.12) 
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and  its  first  derivative 


(4.X3) 


dp(v)  _  y'  dTn(y) 
dy  \  "  dy 

be  evaluated  at  each  collocation  p'  int.  Inserting  p  and  the  other  variables  into  the  linearized 
stability  equations  leads  to  a  system  of  linear  equations  for  the  coefficients  of  the  basis  func¬ 
tions.  From  these,  one  can  either  calculate  a  global  spectrum  using  routines  from  the  IMSL 
library  or  obtain  a  single  eigenvalue  at  a  sequence  of  parameter  values  using  a  continuation 
technique.  Details  are  given  in  Heruorf,  (1990). 

Figures  1-5  all  correspond  to  M, &  ~  2.  The  dependence  of  the  eigenvalues  of  inviscid 
character  is  shown  in  Figures  1  a-c  which  correspond  to  a  =  0.05,0.10  and  0.20  respectively, 
all  at  Re  =  2  x  10s,  and  obtained  with  100  collocation  points.  It  is  clear  from  the  sequence  of 
figures  that  one  inviscid  eigenvalue  is  approximately  stationary,  while  the  remaining  inviscid 
eigenvalues  originate  from  very  large  values  of  Rl{c}  at  very  low  a.  A  magnification  of  figure 
lc  (shown  in  figure  Id)  shows  the  region  of  the  spectrum  corresponding  to  0  <  Rl{c}  <  1. 
One  notes  tne  general  “Y”  shape  of  the  spectrum  as  discussed  by  Morawetz  (1952).  The 
spectrum  is  composed  of  three  “Y”  shaped  pieces.  These  are  due  to  the  structure  of  the 
viscous  terms  of  the  two  momentum  equations  and  the  temperature  equation.  The  coupling 
between  the  equations  and  the  non  unit  Prandtl  number  leads  to  the  non  superposition  of  the 
“Y”  singular  curves  as  Re  -+  oo.  Before  drawing  any  general  conclusions  however,  we  show 
spectra  at  Re  =  7  x  10s  in  Figure  2,  leaving  all  other  parameters  unchanged  with  respect  to 
Figure  1.  As  expected,  higher  resolution  is  required  to  properly  resolve  the  spectrum.  The 
sequence  of  Figures  2a-d  correspond  to  resolutions  of  75,  100,  125  and  150  collocation  points 
respectively.  Although  Morawetz  (1952)  predicts  that  there  is  a  set  of  “viscous”  eigenvalues 
which  lie  on  the  edges  of  the  “V”  part  of  the  “Y”,  packed  with  a  density  proportional  to 
VaRe,  it  is  clear  that  at  M <*,  =  2,  the  triple  point  of  the  “Y”  is  cut  off  by  two  horizontal 
bands  of  eigenvalues.  The  word  horizontal  is  used  here  as  a  qualitative  description.  At 
the  lower  resolution  of  N  =  75,  the  bands  are  slightly  wider  then  for  higher  N.  It  is  not 
completely  clear  whether  these  bands  will  remain  or  disc  ppear  as  N  —*  oo.  A  definitive 
answer  would  require  further  investigation  into  the  properties  of  the  viscous  component  of 
the  spectrum  at  high  Reynods  numbers  and  high  Mach  numbers.  Comparing  Figures  2b  and 
2d,  it  is  clear  that  one  effect  of  insufficient  resolution  is  the  splitting  of  the  spectrum  from 
the  vertical  part  of  the  “Y”  (which  is  indicated  by  “|”).  However,  such  a  splitting  is  not 
always  indicative  of  loss  of  resolution  (see  Figure  5  at  a  =  0.1,  M <*,  =  5,  Re  =  2  x  106).  At 
a  —  3.5,  Re  =  7  x  10s,  the  horizontal  bands  become  more  evident  (Figure  3).  Note  that  the 
density  of  eigenvalues  along  this  band  has  increased  significantly.  Figure  4  shows  a  blow  up 
of  the  upper  left  hand  corner  of  the  stable  diagram.  The  scales  along  Rl{c}  and  Im{c }  are 
identical  to  bring  out  the  real  shape  of  one  of  the  “triple  points”.  At  N  =  150,  one  notices 
the  left  branch  of  the  “Y”  which  is  at  the  theoretical  angle  of  60°.  As  N  changes  from  125 
to  150  nodes,  the  horizontal  bands  shift  slightly  in  the  direction  of  increased  stability.  Thus, 
the  final  converged  state  of  these  bands  is  inconclusive  for  the  chosen  set  of  parameters. 

We  now  turn  our  attention  to  a  similar  set  of  spectra  at  M ©o  =  5  which  seem  to  exhibit 
a  more  complex  character  than  their  M <*>  =  2  counterparts.  At  Re  =  2  x  106,  a  =  0.1,  the 
three  “Y”  shaped  curves  are  clearly  present.  These  structures  subsist  in  an  unmodified  form 
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(i.e.,  without  the  horizontal  bands)  at  higher  Reynolds  numbers  at  the  higher  M One 
“Y”  is  approximately  symmetrical  about  the  Rl{c}  =  0.5  axis.  The  “|”  structure  has  split 
into  three  pieces,  one  along  Rl{c}  =  0.5,  and  the  other  two  curves  placed  symmetrically 
about  Rl{c }  =  0.5.  The  actual  locus  of  points  along  these  curves  shows  a  similarity  with 
the  continuous  temporal  boundary-layer  spectra  discussed  in  Ashpis  and  Erlebacher  (1989). 
As  the  Reynolds  number  is  increased  to  5  x  106,  the  resolution  requirements  simultaneously 
increase.  Only  when  the  vertical  locus  of  eigenvalues  along  Rl{c}  =  0.5  appears,  is  the 
resolution  adequate,  at  least  for  the  part  of  the  spectra  above  it,  although  this  statement  is 
not  quite  true  near  the  horizontal  bands.  The  unchanging  vertical  position  of  these  bands 
as  N  increases  from  125  to  150  indicates  that  its  presence  is  not  an  artifact  of  a  loss  of 
resolution. 

To  complete  the  picture,  Figures  7  and  8  show  the  spectra  at  Moo  =  5,  a  =  3.5  at 
Re  =  2  x  105  and  Re  =  7  x  10s.  In  both  cases  some  inviscid  modes  (those  near  the  Im{c} 
axis)  have  moved  into  the  0  <  Rl{c}  <  1  range.  As  each  mode  crosses  Rl{c}  =  0  or 
Rl{c}  =  1  regions,  a  critical  layer  develops;  this  is  examined  in  Section  6.  Once  again,  as 
the  Reynolds  number  increases,  the  width  of  the  horizontal  bands  increases,  and  resolution 
studies  indicate  that  they  do  not  disappear  as  N  -*  oo. 

In  anticipation  of  the  results  presented  in  the  following  sections,  we  plot  the  evolution 
of  the  mode  I  phase  velocity  (finite  i?/{c}  as  a  -*  0)  as  a  function  of  a.  This  is  carried 
out  for  several  resolutions  and  at  two  different  Reynolds  numbers.  Although  inviscid  theory 
suggests  stronger  instabilities  at  M  =  5.  viscous  calculations  indicate  that  the  resolution 
requirements  become  much  more  severe.  We  therefore  restrict  ourselves  to  M <*>  =  2.  Both 
at  Re  =  7  x  105  and  Re  =  1.4  x  10^.  7m{c}  exhibits  a  local  maximum  near  a  =  4.  This 
mode  becomes  less  stable  at  higher  Re,  while  a  shifts  to  a  slightly  lower  value.  As  expected, 
a  higher  resolution  is  required  at  Re  =  1.4  x  106  ( N  =  150),  as  opposed  to  N  =  100  at 
Re  =  7  x  10s.  Unfortunately,  we  are  not  able  to  increase  Re  beyond  1.4  x  106  and  still  obtain 
accurate  results  at  resolutions  not  exceeding  N  =  150.  Figures  9  and  10  show  the  variation 
of  Rki  and  Rlcr  for  both  Reynolds  numbers.  The  phase  velocity  has  been  scaled  by  \/Re  to 
bring  out  the  one-half  power  scaling  law.  Further  discussion  of  Figure  9  is  deferred  until  the 
asymptotic  theory  at  large  Re  has  been  presented. 

In  the  following  several  sections,  we  examine  solutions  of  the  stability  problem  in  the 
limit  Re  — >  oo.  Guided  by  a  number  of  previous  viscous  theories,  together  with  the  work  of 
Morawetz  (1952,  1959),  we  expect,  in  general,  the  solution  to  develop  in  the  form 

u  =  uo{y)  +  Re~*v,i{y)  4-  0(Re-1), 
v  =  v0(y)  +  Re"^ui(y)  +  0(Re-1), 

f  =  T0(y)  +  Re~^Ti(y)  +  <3(Re-1),  (4.14) 

p  =  P0(y)  +  R^Pi(y)  +  O(Re-1), 
c  ~  cq  +  Re"2ci  +  0(Re-1). 
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5.  Inviscid  Disturbance  Equations 


Taking  equations  (4.3)  -  (4.7),  and  then  retaining  lowest  order  terms  in  Re  (see  Mack 
1984,  for  example),  we  obtain  the  following  two  first-order  equations  for  the  leading  order 
normal  velocity  ( v  =  u0)  and  pressure  (p  =  po)  perturbations  respectively 

-  Uoyvo  ipo  {To  -  MI(Uq  -  cq)2} 


Uo -Co  7 


Uq  —  Cq 


POy  tQ  fTT  \~ 

^r  =  -jr(%-coK 

where  the  appropriate  boundary  conditions  to  be  applied  to  this  system  are 


implying 


uo(0)  =  v0(l)  =  0, 


Poy(0)  =po„(l)  =  0, 


where  Co  is  the  complex  wave-speed  of  this  system.  Equations  (5.1)  and  (5.2)  may  be 
combined  to  yield  the  single  second  order  equation  (see  for  example  Mack  1984,  1987) 


f  (t/p  —  Cq)Vqv  —  ^OyUpl  2 
\  T„  -  -  <*)»  /  " 


(Uo  -  c) 


for  uo- 

Before  investigating  any  of  the  above  systems  numerically,  it  is  interesting  to  study 
the  significance  of  so-called  “generalized  inflection  points”  which  are  found  to  be  highly 
significant  in  the  context  of  compressible  flows.  If  we  multiply  (5.5)  by  Vq  (the  complex 
conjugate  of  Uo)  and  then  subtract  from  the  resulting  equation  its  complex  conjugate,  the 
following  equation  is  obtained: 


up  d  f  (U0 
Jo-c^dyX 


Co)uQy  ~  UoyVo 


\  =  fro  d  [  (Up  -  c*o)VQy  ~  UoyVo 
J  Uo-  Co  dy\  X* 


where 


x  =  To  -  Ml(U0  -  Co)2,  (5.7) 

and  a  superscript  asterisk  here  denotes  a  complex  conjugate.  After  some  algebra,  this  may 


be  written 


d  fuoyl  _  d 


u0j 

dy  [  x 


Voy  _  -  f_l  d  Upy _ 1  d  Upy  )  /K  o\ 

X*J  ~  V°Vo\U0-cody  x  Uq  —  c^dy  [  X*  J  J  ' 


Writing 


co  =  /?/{co}  -f  i  Jm{co} 

then  the  neutral  state  corresponds  to  7m {cq}  — >  0.  In  this  limit,  (5.8)  may  be  written 


(5.9) 


(5.10) 


'VqVOv  -  V0Vov 

vi  Vov(yi)2ri  / 

x 

Vo  UoiiiVi)  ' 

U0y\ 

T^- 


Using  arguments  similar  to  those  of  Duck  (1990),  then  (i)  as  /m{co}  -+  0,  the  left-hand-side 
of  the  equation  is  zero  except  possibly  at  (where  U0(yi)  ==  /?«’{co}),  (ii)  the  term  in  the 
parentheses  on  the  left-hand-side  of  (5.10)  must  be  zero  at  both  y  =  0  and  y  =  1,  while  (iii) 
the  right-hand-side  is  clearly  non-zero  unless 


£  m 

dy  To  J  y-v. 


=  0. 


(5.11) 


Thus,  in  order  to  avoid  an  inconsistency,  we  must  have  (5.11). 

Most  importantly,  there  does  exist  a  difference  between  the  present  (bounded)  flow  con¬ 
figuration,  and  that  of  unbounded  flows,  in  that  in  the  present  situation,  (5.11)  is  necessary 
if 

0  <  Rl{co)  <  1,  (5.12) 

i.e.,  a  critical  layer  must  exist  inside  the  flow,  while  in  the  case  of  unbounded  flows  condition 
(5.11)  only  holds  if  the  wavespeed  is  “subsonic”,  i.e.,  1  -  1/M^  <  Rl{co}  <  1  +  1  /Mw.  This 
is  because  (5.11)  is  a  direct  consequence  of  the  zero  velocity  perturbations  at  the  domain 
boundaries.  There  is  absolutely  no  distinction  made  here  between  supersonic  and  subsonic 
modes.  Thus  any  neutral  inviscid  mode  satisfying  (5.12)  must  be  associated  with  a  general¬ 
ized  inflection  point.  The  above  says  nothing  about  neutral  inviscid  modes  outside  the  range 
of  (5.12).  It  is  also  worth  noting  that  authors  who  implement  Dirichlet  boundary  conditions 
(in  place  of  radiation  boundary  conditions)  on  the  disturbance  terms  in  truncated  infinite 
domains,  may  well  experience  difficulties  computing  non-inflectional  supersonic  modes,  since 
the  arguments  above  suggest  that  a  generalized  inflection  point  is  necessary  for  supersonic 
disturbances,  a  condition  which  is  clearly  erroneous  in  the  unbounded  case. 

Utilizing  (3.5)  and  (3.9)  in  (5.11)  yields, 


=  0.  (5.13) 

If  the  term  inside  the  square  brackets  is  zero,  then  by  Sutherland’s  law  (2.9),  we  must 
have 

[37b  +  5C]y=y.  =  0,  (5.14) 

which  is  clearly  inadmissible.  Consequently,  the  only  way  that  a  generalized  inflection  point 
will  occur  is  if  the  mean  temperature  profile  has  a  local  extremum.  If  we  invoke  Sutherland’s 
law  and  (3.9),  then  we  either  require 


UoyiVi)  —  0) 

(5.15) 

which  is  clearly  not  possible  on  account  of  (3.5),  or 

rj  („  x  _  (1  ~  r)Tr 

Uo{yt)  2pnr  _  1]  • 

(5.16) 
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Since  Tr>  1,  and  the  flow  is  unidirectional,  then  this  condition  cannot  be  satisfied  unless 


which  implies 


0  <  Uo(yi)  <  1, 


1  _  2=±oM2 
- - <r<l, 


1  + 

(implying  that  the  lower  wall  must  be  cooled  below  adiabatic  conditions). 


(5.17) 

(5.18) 


6.  Inviscid  Disturbance  Results 

For  this  set  of  computations,  the  basic  profile  for  velocity  and  temperature,  namely  Uo(y) 
and  To(y)  was  determined  by  means  of  a  fourth-order  Runge-Kutta  scheme.  The  eigenvalue 
problem  posed  in  (5.1)  -  (5.4)  was  solved  similarly,  with  Newton  iteration  being  used  to 
update  the  complex  wavespeed  Co  so  that  all  boundary  conditions  were  satisfied. 

In  a  number  of  computations  (specifically  those  for  which  0  <  Rl{co}  <  1),  it  was  found 
necessary  to  deform  the  integration  contour  into  the  complex  y  plane  (in  particular  below 
the  real  y  axis),  in  a  manner  described  by  Lees  and  Lin  (1946).  This  was  undertaken  by  two 
independent  ways.  The  first,  as  suggested  by  Mack  (1965b),  involves  obtaining  the  mean  flow 
profile  (Uo(y),To(y))  along  the  real  y  axis,  and  then  using  Taylor  series  expansions  to  obtain 
the  mean  flow  profile  off  the  real  y  axis;  this  detour  is  made  close  to  where  U0(y)  =  Cq.  The 
second  approach  (which  was  generally  used  in  preference)  obtains  the  mean  flow  solution 
itself  in  the  complex  y  plane,  thereby  eliminating  the  errors  associated  with  truncation  of 
the  Taylor  series.  Specifically,  the  mean  flow  was  obtained  for  0  <  y  <  yi  (i.e.,  on  the  real  y 
axis),  yi  <  y  <  yi  -  iy 2,  yi  -  iy?  <y  <y$  —  iyj,  Vs  -  iy 2  <  y  <y  3  and  thereafter  back  along 
the  real  axis  y3  <  y  <  1.  yi,  2/2?  2/3  were  all  taken  to  be  real  and  positive  (although  this  is 
not  essential)  and  were  chosen  to  avoid  the  computation  proceeding  too  close  to  the  critical 
layer.  Comparison  of  results  using  the  two  approaches  proved  a  useful  check  of  the  accuracy 
of  our  results.  A  further  check  on  our  results  was  that  in  addition  to  solving  (5.1),  (5.2)  we 
also  solved  the  adjoint  system  (see  (7.29),  (7.30)  below),  and  also  (5.5). 

The  first  results  we  present  is  for  the  case  M, =  2  and  adiabatic  lower  wall  boundary 
conditions  (here  and  in  all  cases  we  took  7  =  1.4,  a  =  0.72  together  with  the  Sutherland 
constant  C  =  0.5).  Results  for  RI{cq}  are  shown  in  Figure  11.  The  results  display  a  number 
of  interesting  features.  There  appear  to  be  many  modes,  which,  however,  can  be  divided 
into  two  distinct  families,  the  first  corresponding  to  Rl{co }  >  1  as  a  — ►  0.  All  these  modes, 
with  the  exception  of  one  (mode  I)  have  i?/{co}  — >  00  as  a  — >  0;  we  refer  to  these  as  the 
upper  family.  These  modes  all  suffer  a  monotonic  decrease  in  RI{cq}  as  a  increases,  and 
the  results  suggest  that  some  finite  value  is  approached  in  the  limit  a  —*  00.  The  second 
family  is  defined  by  Rl{co}  <  0  as  a  — ►  0,  and  the  Rl{co}  all  increase  monotonically  with  an 
increase  in  a,  and  eventually  become  positive.  We  refer  to  these  modes  as  the  lower  family 
of  solutions.  Again,  all  except  the  first  of  these  modes  (mode  II)  appear  to  be  unbounded 
as  a  — »  0. 

A  further,  important  feature  is  that  in  all  cases  for  which  Rl{co)  >  1  or  i?/{co}  <  0,  all 
these  modes  are  neutrally  stable  to  within  machine  accuracy,  i.e.,  Im{co}  =  0.  However, 
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since  there  can  be  no  generalized  inflection  point  for  the  mean  flow  under  consideration,  then 
there  can  be  no  neutral  modes  with  0  <  cq  <  1.  Indeed,  it  is  found  that  in  the  case  of  the 
upper  family  of  solutions,  once  7?/{co}  <  1,  then  these  modes  become  stable  according  to 
our  inviscid  calculations,  while  in  the  case  of  the  lower  family  of  solutions,  once  Rl{co}  >  0, 
then  these  modes  become  unstable  according  to  our  inviscid  calculations;  the  location  at 
which  Rl{co}  either  drops  below  1  or  rises  above  0  is  marked  on  Figure  11  by  a  circle.  The 
distribution  of  7m{co}  for  mode  I  is  shown  in  Figure  12  (other  upper  family  modes  also 
have  negative  values  of  7m{co),  but  many  orders  of  magnitude  smaller).  The  distribution  of 
7m{co}  for  mode  II  is  shown  in  Figure  13  (other  modes  of  the  lower  family  have  considerably 
smaller  values  of  7m{co},  although  these  are  also  positive).  Next  we  go  on  to  consider  the 
aforementioned  trends  as  a  — >  0  suggested  by  our  numerical  results. 

The  existence  and  behavior  of  modes  I  and  II  as  a  — ►  0  is  easy  to  confirm.  If  we  set 
a  —  0  in  (5.5),  and  integrate  once,  we  obtain 


(Uo  -  cc)v0y  -  UoyV o  =  K[T0  -  Mi(U0  -  Co)2],  (6.1) 


where  K  is  an  arbitrary  constant.  This  equation  may  be  integrated  once  more  to  yield 


vo 


I<(Uo-co)  J“ 


{T0  -  Ml(U0  -  an 

- (£^p - V 


(6.2) 


If  this  is  to  satisfy  the  boundary  condition  on  y  =  1  (the  boundary  condition  on  y  =  0  is 
already  satisfied  by  (6.2)),  we  must  then  have 


£(U^rdy=Ml  (63) 

This  relationship  was  solved  for  Co  (for  given  Uo(y),  T0(y),  and  M^)  and  confirmed  the 
numerical  results  as  a  — >  0.  Indeed,  if  we  take  a  model  problem,  i.e., 


Uo(y)  =  y,  T0(y)  =  1  (6.4) 

then  the  solutions  of  (6.3)  are 

.  1  ±  [1  +  4 /Ml]'/’ 

a  = - 2 - •  (6-5) 

These  modes  will  be  discussed  further  in  Section  7. 


The  behavior  of  modes  III  and  higher  as  a  —*  0  is  different,  but  nonetheless  straightfor¬ 
ward  to  confirm.  In  these  cases  we  have  co  — >  oo  as  a  — »  0,  and  a  balancing  of  terms  in  (5.5) 
in  this  limit  demands 


c0  =  a-1  coo  +  0(1), 

(6.6) 

vq  =  5oo  +  0(a). 

(6.7) 

The  equation  for  vqo  is  then 

~  .  cIqM^voo  n 

Voyy  +  =  u- 

(6.8) 
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If,  to  make  further  progress  analytically  we  assume  |coo|  1,  implying  we  are  focusing 
attention  on  the  higher  nodes  (alternatively  we  could  assume  the  model  problem,  with  To  = 
1),  then  the  WKB  solution  is 

Voo  =  AtJ  je<CooM“  fo  T^y  -  e-<cooMoo  /0y r0"1/a^|  ?  (6.9) 

where  A  is  a  constant.  Here  the  boundary  condition  on  y  =  0  has  been  applied,  while  the 
other  boundary  condition  on  y  =  1  demands 


CqO  = 


rnr 


(6.10) 


where  n  is  any  (large)  positive  or  negative  integer.  This  expression  clearly  illustrates  the 
multiplicity  of  modes. 

The  second  set  of  results  we  show  is  for  the  higher  Mach  number  case  M*,  =  5,  with 
the  adiabatic  boundary  conditions  on  the  lower  wall  retained.  Figure  14  shows  a  number  of 
results  for  Rl{co}.  Although  these  seem  much  the  same  as  the  corresponding  distributions 
for  Moo  =  2,  there  are  some  differences,  in  particular  with  mode  I  a,  as  defined  in  Figure 
14,  which  initially  corresponds  to  mode  I  in  Figure  12.  However,  unlike  the  corresponding 
results  for  Moo  =  2,  it  turns  out  that  there  exist  other  (stable)  modes  in  the  vicinity  of  mode 
1a-  Some  of  these  modes  are  shown  in  Figure  15  (on  a  magnified  scale)  and  are  denoted  by 
Ib  and  Ic-  These  modes  were  difficult  to  generate  numerically  in  regions  of  a  where  RI{cq} 
was  very  close  to  unity  and  where  the  decay  rates  were  quite  small;  however,  the  existence 
of  these  modes  was  confirmed  using  the  three  different  formulations  of  the  inviscid  stability 
problem  described  above,  together  with  the  condition  (6.3).  Indeed,  it  is  quite  likely  that 
other  such  modes  exist,  although  the  present  investigation  failed  to  yield  any.  However,  it 
turns  out  that  these  modes  are  all  stable,  as  shown  in  Figure  16.  Note  that  modes  1a  and 
Ic  have,  in  places,  values  of  cq  exceedingly  close  to  each  other,  but  nonetheless  distinct;  this 
compounded  the  difficulty  associated  with  these  modes. 

Modes  II,  and  higher,  on  the  other  hand  did  exhibit  the  same  qualitative  behavior  as  the 
corresponding  M0 0  =  2  results.  In  particular,  mode  II  becomes  unstable  at  a  ~  1.85,  with 
the  distribution  of  Im{co}  shown  in  Figure  17.  Other  higher  modes  of  this  lower  family  are 
also  unstable,  but  with  substantially  smaller  growth  rates. 

The  final  set  of  results  in  this  section  relates  to  the  case  M^  =  5,  as  previously,  but  with 
the  (cooled-wall)  boundary  condition  7o(0)  =  1.  It  is  worth  noting  that  even  in  this  case,  the 
uniform  temperature  together  with  a  linear  velocity  profile  is  still  not  a  proper  solution  of 
the  governing  equations,  except  in  the  very  special  limit  as  7  — >  1.  Nonetheless,  the  solution 
does  have  some  important  properties,  namely  that 


T0(y)  =  T0(l  -  y),  (6.11) 

together  with 

U0(y )  =  1  -  Uq(1  -  y).  (6.12) 

As  a  consequence  of  these  symmetries,  the  mean  flow  has  a  single  G.I.P.  at  y  =  §,  (with 
Uq  =  |).  These  properties  turn  out  to  have  interesting  implications  for  the  inviscid 
stability  of  the  profile. 
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Figure  18  shows  the  variations  of  Rl{co)  for  the  first  eight  modes.  As  in  the  previous  cases, 
there  again  exist  two  distinct  families,  with  all  except  two  of  the  modes  having  the  property 
\Rl{co}\  — *  oo  as  a  — »  0.  Furthermore,  all  modes  are  neutrally  stable  if  i?/{co}  <  0  or 
Rl{co}  >  1.  Significantly,  due  to  the  aforementioned  symmetries  about  y  =  |,  the  behavior 
of  Rl{co}  for  the  upper  family  of  modes  is  merely  the  mirror  image  of  the  corresponding 
lower  family  member  about  y  =  These  symmetries  also  yield  the  result  that  Im{c 0}  is 
precisely  the  same  for  corresponding  modes  at  the  same  wavenumber.  Figure  19  shows  the 
distribution  of  7m{co}  for  modes  I  and  II.  These  modes  begin  as  neutrally  stable,  then  as 
Rl{co}  drops  below  unity/rises  above  zero,  both  modes  become  stable.  As  a  increases  still 
further,  both  modes  have  RI{cq}  —►  |  which  is  reached  at  a  «  5.5.  Jm{co}  then  becomes 
positive,  implying  unstable  modes.  This  is  entirely  consistent  with  our  G.I.P.  arguments, 
which  predict  neutral  modes  having  Co  =  5.  Thereafter,  RI{cq)  for  mode  I  increases,  and 
decreases  for  mode  II. 

Figure  20  shows  the  distribution  of  /m{co}  for  modes  III  and  IV,  which  over  the  range 
shown  are  seen  to  be  stable.  It  is  possible  that  regions  of  instability  exist  for  higher  a, 
although  the  computations  become  increasingly  difficult  as  a  — *  oo. 

In  the  following  section,  we  go  on  to  consider  regions  where  viscous  effects  are  likely  to 
become  important,  particularly  in  determining  the  stability  properties  of  the  flow. 


7.  The  Viscous  Correction 


The  results  described  in  the  previous  sections  show  that  there  is  a  family  of  solutions  of 
(5.1)  -  (5.4)  comprising  neutrally  stable  modes  over  a  wide  range  of  a,  which  have  either 
Co  >  1,  or  c0  <  0,  and  hence  have  no  critical  layer.  The  question  then  arises  as  to  the  effect  of 
viscosity  on  these  modes  -  whether  it  plays  a  stabilizing  or  destabilizing  role  in  this  problem. 
We  investigate  this  question  next. 

The  viscous  correction  to  the  problem  arises  due  to  the  thin  layers  that  occur  as  a  result 
of  the  violation  of  the  no-slip  condition  (and  also  the  temperature  condition)  on  the  walls 
y  =  0  and  y  =  1.  These  conditions  are 

tt  =  t>  =  0  on  y  =  0  and  y  =  1,  (7.1) 


T  =  0  on  y  =  1,  while  on  y  =  0, 

dT  n 

o-  =° 

oy 

(7.2) 

if  the  wall  is  insulated,  and 

o 

II 

(7.3) 

if  it  is  heated  or  cooled. 

Now  from  (4.4)  and  (4.6)  we  see 

that  on  y  =  0  we  have 

...  To(0)po(0) 

yM^co 

(7.4) 

and 

=  r„(o)p(o), 

To 

(7.5) 

15 


while  on  y  =  1  we  have 


r,(i)fe(i) 

0  7MJ,(cq-1)’ 


(7.6) 


and 

To  =  Io(l)po(l).  (7.7) 

Thus,  generally  these  expressions  fail  to  satisfy  the  appropriate  wall  conditions  described 
above. 

In  order  to  remedy  this,  we  require  a  (Stokes-like)  layer  of  thickness  0(Re-'*)  on  both 
walls.  Considering  first  the  layer  on  y  =  0,  defining 


Y  =  t/Re*  =  0(1), 

(7.8) 

togethef  with 

u  =  U(Y)  +  0(  Re~*), 

v  =  Re~*V(Y)  +  0(Re“1), 

f  =  f{Y)  +  0{  Re"*), 

(7.9) 

p  =  P(Y)  +  0(  Re"*), 

p  =  p(F)  +  0(Re“*), 

then  taking  0(1)  terms  in  (4.3)  -  (4.7),  we  have 

—iacoU  ~  iaP 

X„(0)  ~>‘lVyy  tm*,’ 

(7.10) 

II 

o 

(7.11) 

•  X  ,  iu  ,  1  f, 

,Co'’+r„(o)  +  r(o)F,'  =  0’ 

(7.12) 

-  ^Trr  +  (— )  iacoPo  =  0, 
To(0)  C7  \  7  / 

(7.13) 

(7.14) 

where  pi  =  po{0).  On  account  of  (7.11)  the  pressure  must  be  constant 
so  we  write 

across  this  layer,  and 

P  =  Po(°)  =  Pl- 

(7.15) 

The  solution  to  (7.10)  is  then 

0_PlT b(0)  f,  ,-(,»)  M 

(7.16) 

7M£>co 
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In  the  case  of  an  insulated  wall  on  y  =  0,  the  temperature  perturbation  is  given  by  (7.5) 
across  the  layer  (i.e.,  it  is  constant),  while  for  the  heated  or  cooled  wall  conditions,  using 
(7.13)  we  obtain 

m  _■( 7  ~  ~  rrifn\  i  1  ~  (iir  Tn?0))  ^  1  (n  1 


-)piT-(O)  |l  - 


(7-17) 


Using  the  continuity  equation,  (7.12),  we  obtain 


V  =  ico PlY  -  icopL 


y  + - - - .-ISSft)  V-1 

—iacnff  1  2 

[mlTo(0)J  L 


(7.13) 


.!Mo(0)Jy+ - I -  c  (’Smfe)  Y  -1  l, 

IfM^Co  I  f  -tea  1  a  | 


where  6  =  0  in  the  case  of  an  insulated  wall,  while  b  =  1  for  a  cooled/heated  wall.  The 
corresponding  layer  on  the  upper  wall  is  quite  similar.  Defining  the  lengthscale 

y  =  (S  -  l)Re*. 

together  with  (7.9),  then  the  results  for  Y  =  0(1)  may  be  simply  inferred  from  those  for 
Y  =  0(1)  by  (i)  replacing  ‘F’  by  Y\  (ii)  replacing  ‘pi,’  by  lpu'  =  p(l),  (iii)  replacing  lCo’ 
by  “co  —  1”  and  (iv)  replacing  \il  and  jTo(0)  by  unity.  Finally  the  (key)  result  for  V  is  then 
given  by 

V  =  i(co  -  1  )pvY 


-i(co  -  1  )pu 


7-1 


Y _ _ _ _  [eH°(co-iW*v  _  j 

[-m(co  -  1  )cr]a  *■ 


(7.19) 


7  Ml(co-l) 


-ia(co  -  1)] 


- -  fe[-Mco-i)]^  _  i]  \ 

-111*  L 


We  shall  later  require  the  behavior  of  (7.18)  and  (7.19)  as  7  -+  oo  and  Y 
respectively.  We  have 

V  Y  l icopL  -  icoPL  f ^ 


,  icoPLb( _  ipLTo{0) 

'  r  •  +  v 

-iaco<r  2  'vM2  Cn  ~tac0-  ■ 
To(0)w  7Mco°0  To(0)hl. 


(7.20) 
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vqlY  +  vil  as  Y  -*  oo, 


(7.21) 


V  -  yUc- 

l  f-«a(co- 


[-t£»(co-l)ff)i  tM»(co— 1){— »'«(co— l)]s 


»(co-l)pt/(^Ij _ _ 

[-•a(eo-l)(r]i/4'  7Mj,(eo-l)[-ta(co-l)l^4 

A  A 

— >  tiof/K  +  Uic/  as  y  — +  —oo. 


(7.22) 


(7.23) 


We  now  seek  the  leading-order  viscous  correction  term  to  (5.1),  (5.2)  for  y  —  0(1), 
using  the  expansions  (4.14).  The  variables  iii  and  pi  turn  out  to  be  determined  by  an 
inhomogeneous  form  of  (5.1)  and  (5.2),  namely 


*  EVi  »Pi  \Tq-MKUo-cq )2] 


Uq  —  Cq  Uq  —  Cq 


ia27M4  ...  .  _ 

Ply  +  - ~ - (Uq  —  Co)Ui  =  C1R2, 

■lo 


—  c\R\ 


where 


n  Uo»*°  1  ipoTg _ ipo 

1SS  (Uo-co)>  +  iM^Uo-co)**  7’ 


i'yM^c^vo 


(7.24) 


(7.25) 


(7.26) 


(7.27) 


The  boundary  conditions  on  this  system  are  also  inhomogeneous,  namely 

=  !“>  (7.28) 

Ul(l)  =  Viv, 

(the  boundary  conditions  for  piy  may  be  obtained  from  (7.25),  although  this  is  not  necessary 
for  the  following),  where  v\i  and  vw  are  defined  by  (17.21)  and  (7.23)  respectively. 

To  obtain  ci,  we  use  the  condition  of  solvability  of  the  system  (7.24)  -  (7.28).  For  this 
we  require  the  adjoint  to  the  system;  if  we  denote  v+  and  p+  as  the  adjoint  functions,  then 
these  are  to  be  determined  by 


CV+  «.»((/,  -  cobAfe 
"»+t/„-co  +  To  P 

+  i v*  [To-  Ml(Vo-co)2]  „ 

p’  ~  Uo-co  L — Wi — J“°' 

The  boundary  conditions  to  be  applied  to  this  system  are  that 

p+(0)  =p+(l)  =  0. 


Ci  is  then  given  by 


__  u+(l)uu/  -  u+(0)i;iL 

“  £  [/?iu+  +  R2p+}  dy  ’ 


(7.29) 


(7.30) 


(7.31) 


(7.32) 


where  v*  and  p*  are  the  complex  conjugates  to  v+  and  p+  respectively. 

Given  the  nature  of  the  solution  for  u+,p+,t'it/,uiL,  i?,  and  R2  when  cq  is  real,  it  is  easy 
to  show  that  we  must  have 

ci  =  /m{Cl}(±l  +  »),  (7.33) 

where  the  positive  sign  is  taken  for  cq  >  0,  the  negative  sign  is  taken  for  Co  <  0;  note 
that  /m{ci)  >  0  for  instability.  However,  in  all  the  computations  we  performed,  without 
exception,  we  found  Im{ci}  <  0.  Results  for  7m{ci)  for  the  adiabatic  case  with  M, =  2 
(modes  I  and  II)  are  shown  in  Figure  21.  The  distribution  becomes  unbounded  as  a  — *  0, 
and  also  as  a  —*  era,  the  point  at  which  Co  — ►  1  or  0;  these  distributions  are  typical  of 
all  modes.  Figure  22  shows  the  corresponding  distributions  for  the  adiabatic  case  with 
A7oo  =  5,  and  exhibit,  quai'tatively,  the  same  behavior.  (The  corresponding  computation  for 
M-o  =  5,  Tq( 0)  =  1  was  also  carried  out,  and  yielded  qualitatively  the  same  picture.)  We 
shall  defer  discussion  of  the  limit  a  — ►  a0  until  the  following  section,  but  let  us  now  consider 
the  limit  as  a  — *•  0,  of  modes  I  and  II  from  which  it  appears  that  our  expansions  cease  to  be 
uniformly  valid.  This  is  most  clearly  seen  by  the  apparently  singular  behavior  of  the  viscous 
correction  to  the  complex  wavespeed  as  a  -*  0,  together  with  the  0(a~a)  growth  of  the 
wall-layer  thickness  in  this  limit  (see  (7.10)  for  example).  Thus  when  a  =  0(Re_1)  the  wall 
layers  will  fill  the  entire  channel,  and  the  disturbances  become  purely  viscous  in  nature. 

Specifically,  if  we  write 


p 

II 

% 

CD 

1 

M 

PI 

(7.34) 

then  the  governing  equations  (4.3)  -  (4.6)  reduce  to 

—icp  +  ipoii  +  pov  +  iUop  +  vpoy  =  0, 

(7.35) 

r  .  „  .  ...  „  „TT  ,  tap 

ctpo\—icu  +  iUqu  -f  vUoy  +  —rpr 

7M£ 

(7.36) 

Pottyy  4"  Po-rTUoy  4*  ^orTbyity, 

Py  =  0, 

(7.37) 

apo[~icf  -f  iUoT  +  uToy]  —  [iaUo  —  iac]p 

, rp  |  f^OT  rp  ;ji  ,  ^ f^OT  rp  ,  ^yf^QT 

„  lyy  '  _  1o yJ-y-r  J-o yy -r 

<7  £7  <7  <7 

(7.38) 

+T  /j.oTTT'oy  +  Tp  or(7  —  l)M^,I/oy, 

where  it  has  been  implicitly  assumed  that  viscosity  is  a  function  of  temperature  only.  The 
problem  then  reduces  to  a  basically  viscous  system. 

However,  for  modes  III  and  higher,  cq  =  0(a~1)  as  a  — >  0,  and  since  the  Stokes  layer 
thickness  remains  0(R~ »),  it  does  not  fill  the  entire  channel.  Additionally,  as  a  — >•  0, 
7m{ci}  =  0(a_1),  and  so  the  expansions  (4.14)  remain  valid. 

Returning  to  Figure  9,  we  note  the  general  excellent  agreement  between  the  values  of 
y/ReIm{c }  obtained  numerically,  and  7m{ci}  obtained  from  the  above  asymptotic  theory, 
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for  the  case  M <*,  =  2  (mode  II).  The  values  of  c,  obtained  using  the  two  approaches  are 
almost  indistinguishable.  We  do  note  a  deviation  between  the  results  as  Cq  — >  0,  caused 
by  the  (expected)  breakdown  in  the  asymptotic  theory  in  this  limit  (and  as  cq  — ►  1).  This 
aspect  is  investigated  in  the  following  section. 


8.  The  Nature  of  the  Solution  as  Co  ->  0  (or  cq  — ►  1) 

We  consider  here  the  nature  of  the  complex  wavespeed  as  Co  — >  0  (the  results  for  Co  — +  1 
may  be  simply  inferred  from  those  of  cq  — *  0).  Our  previous  results  indicate  that  as  c0  — >  0, 
(with  a  — ►  c*0)  then  (i)  Cj,  the  viscous  correction  to  the  wavespeed  becomes  unbounded, 
(ii)  the  leading  order  wavespeed  Co  becomes  non-neutral  as  a  increases  above  a0,  and  (iii) 
the  thickness  of  the  wall  layer  increases.  Thus  this  region  is  likely  to  be  a  regime  of  some 
interest,  to  which  we  now  focus  our  attention. 

More  specifically,  guided  by  our  numerical  results,  if  we  write 

a  =  <*0  +  £<*i>  (8.1) 


where 


e  =  Re“», 


then  we  expect 

c  =  eci  +  O(e^).  (8.3) 

The  solution  in  the  “core,”  i.e.,  away  from  the  wall  layers  is  then  expected  to  develop  in  the 
form 

v  =  v0  (y)  +  ei>i{y)  +  •  •  • ,  (8.4) 

P  =  Po(y)  +  epi(y)  +  •  •  •  .  (8.5) 

The  leading  order  system  is  then 

UoyVo  ipo _ »y2  rr2l  _  n  /O  c\ 


or  symbolically 


*  \rp  tt2]  a 

0y  U0  'iMlUo^0  MM  ’ 


Ci{vo,po}  =  0, 


or  symbolically 


-hi-+ia  2^1  =  0 

7  To 


^{voiPo}  -  0) 


together  with 

wo(0)  =  w0(l)  =  0.  (8.10) 

Notice  that  (8.6)  and  (8.8)  imply  that  p0  =  0(y3)  as  y  ->  0;  this  system  effectively  determines 
the  value(s)  of  a0  for  which  cq  =  0.  Turning  to  the  next  order  system,  we  find 


£.{«.,  A}  =  Uo,vj^  +  ,Pi  +  MluS)  =  c ,A„ 


(8.11) 
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C^UPl}  = 


ialcivo  2iaQOiiUoVo 


—  C\R2  +  OC1R3. 


The  boundary  conditions  to  be  applied  to  this  system  are 

u1(l)  =  0,  «!(0)  =  A MO), 

where  A  is  to  be  determined  later.  Using  conditions  of  solvability,  we  must  have  that 

*  -vmApm-oifo'&P+dy 


IW+  +  R2p+]dy 

where  p+(y),v+{y )  are  the  adjoint  functions  to  (3.6),  (8.8),  namely  those  determined  by 
(7.29),  (7.30)  with  a  =  ao,  cq  =  0. 

We  now  consider  the  effect  of  the  wall  layers.  The  upper  wall  layer  remains  of  thickness 
0(Re-^)  and  as  such  plays  no  role  to  this  order.  The  interesting  changes  are  related  to  the 
lower  wall  layer,  where  now  we  must  have 


Y  =  e'xy  =  0(1) 


(8.15) 


as  the  crucial  scale,  wherein  to  leading  order 


(uyv,p)  =  {U0,eVo,eP0).  (8.16) 

Taking  the  leading  order  terms  in  the  continuity  and  momentum  equations  we  find 


iUo  +  V0y  =  0, 


(8.17) 


P0Y  =  0,  (8.18) 

iao[/7o(0)y  -  cijUo  +  ao^o(0)Vo  =  U0yy  ~  ia0Po,  (8.19) 

with  boundary  conditions 

£/o(0)  =  Vo(0)  =  0, 

and  as  Y  — >  00. 

Uq  ->  -iv 05,(0),  (8.20) 

where  the  latter  condition  arises  from  a  proper  matching  with  the  y  =  0(1)  solution.  If  we 
now  differentiate  (8.19)  with  respect  to  Y,  invoke  (8.17)  and  (8.18)  we  find 

Uqyyy  ~  *'«o[t/'(0)r  -  cx\UQy  =  0,  (8.21) 

whose  solution  is 

U0y  =  ^AiHiaoU^mHY  -c1)}, 

where  C\  is  a  constant  to  be  determined,  and  cx  =  Ci/Uq(0). 

Utilizing  the  boundary  condition  on  Y  =  0,  and  also  evaluating  (8.19)  on  Y  =  0  leads  to 

iaoPotf  Ai{[mo^(0)]*(y  -  cx)}dY 


[z«o£^(0)]  ^  -  [i«o^(0)  j  3  ca } 


Vo  may  now  be  determined  by  continuity  (8.12)  together  with  (8.20),  and  is  given  by 

p  Po  f?  ft  At{[tao^(0)]i(y  -  c1)}dYdY 
°  “  [iaoU^Ai'i-liaoU^hx) 

Correct  matching  with  the  core  then  requires 

-K£7i(0))i«o,(0)>li'{-[ioo^(0)]bI} 


Po  = 


together  with 
and 


el  Jo“  -4i{[ia„l/i(0)]i(K  -  c,)}<iy 

A =pi(o), 


(8.23) 


[Vi(n  -  Wrf.  (K))^  =  6,(0). 

It  is  these  conditions  that  then  enable  us  to  determine  A(ci)  introduced  previously, 


<*l  C  £  At{[iaoE/j(0)]*(y  -  c,)d?dY 
(ia£/i(0)]Ui'{-[iao£7i(0)]ic,} 


A  = 


(8.24) 

(8.25) 

(8.26) 

(8.27) 


The  procedure  is  then  to  determine  a  solution  to  (8.11),  (8.12)  consistent  with  (8.13),  (8.14), 
and  (8.27).  However,  the  problem  for  cj  is  highly  nonlinear,  and  as  such  is  a  difficult 
numerical  task.  However,  if  we  suppose  |cti|  — ►  oo,  then 


Cl  = 


— fo  Ap+dy 

Jo1[&0+  +  A2p+]dJ, 


0(a"*), 


(8.28) 


which,  on  account  of  the  nature  of  p+,  v+,  i?i,  i?2,  #3  can  be  shown  to  be  purely  real  to 
leading  order:  indeed  it  is  interesting  to  note  that  although  this  sytem  does  admit  complex 
values  of  ci  for  au  =  0(1),  it  appears  that  the  above  fails  to  capture  the  non-neutral  nature 
of  the  inviscid  modes  for  a  >  ao- 

However,  the  reason  for  this  is  clearly  illustrated  by  considering  the  nature  of  the  inviscid 
system  (5.1)  -  (5.2)  in  the  limit  as  a  — >  ao.  Suppose  that  we  set 

a  =  ao  +  a,  (8.29) 


where 

|6|  <  M,  (8-30) 

then  we  expect  expansions  of  the  form 

no  =  Vo  +  olv i  J-  a2v2  +  0(a3), 

Po  =  Po  +  ctpi  +  o?f>2  +  0(d3), 
co  =  aci  +  d2c2  +  +0(a3).  (8.31) 
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The  leading  order  system  is 


A{*o,Po}  =  o, 


(see  (8.7),  (8.9)),  with 


A{vo,p0}  =  0, 


«o(0)  =  60(1)  =  0. 


(8.32) 

(8.33) 


At  the  next  order,  Vi  and  pi  are  again  determined  by  means  of  (8.12),  although  we  shall 
defer  discussion  of  the  boundary  conditions  to  be  applied  to  this  system,  and  turn  instead 
to  consideration  of  the  0(d2)  terms.  We  find 


£i{t)2,p2} 


UoyV^Ci  ,  v0U,vc2  t  UoyV\C\  t  ipiCi  fn,  „r2  rr2l 

~vf  ■ + ~ur  ~w  ~uf[  ° "  M 


ip 0 
+  C/o2 


*  ,  c* 
Cl+u0 


(To  -  MlUl ]  +  2 MlU0c\  +  2Mlcm 


(8.34) 


=  R 


o? 


C2{i)2,p2}  = 


i  a  2 toto  a  ,  2ta0ci  a 

-7frU0vo  -  -=-UoV\  +  —^—v0 
j-  0  -*0  -*0 


i® 0  f  »  »  A  A  I 

—  [-CxVj  -  C2l>oJ 


=  Rb- 


(8.35) 


However,  the  expansions  above  are  not  uniformly  valid  as  y  -*  0,  since  U0(y)  —  0(y) 
in  this  limit,  and  hence  a  breakdown  to  our  approximations  must  occur,  specifically  when 
y  —  0(a).  If  we  set 

y  =  &7,  Y  =  0(1),  (8.36) 

then  Vo  develops  in  the  following  manner  (that  may  be  demonstrated  to  be  correct  a  poste¬ 
riori) 

vo  =  d$0  (F)  +  a2  In  d$x(F)  +  d2$2(F)  +  •  •  •  (8.37) 

It  is  simple  to  show  that 

*o(F)  =  F,  (8.38) 

Ci 

where  K0  is  some  constant;  this  is  consistent  with  (8.33).  We  shall  defer  discussion  of  $i(F) 
until  later,  while  $2(F)  is  given  by 

$2  7 _ t%(0)*2  =  KlUO)  KoYnjO) 

MW  -  Cl)  [0o(O)F  -  exp  [^(0)K  -  C!]2  +  ro(0)[^(0)F  -  cxp 

__  (8.39) 

c2Aq2q(0)  \Uq(0)Y2Ko 
Ci[ WF-cxP  ai[W)y-ci]2’ 

where  K\  is  a  constant. 


23 


Asymptotic  analysis  of  (8.39)  reveals  that  as  F  — ►  oo 


$2  A0F  InF  +  AXF  •+■  A2  In  F  +  A3  +  •  •  • ,  (8.40) 

where  the  A,-  are  all  constants.  It  is  easy  to  show  (see  below)  that 

\  KoTW)  UmKoCi  (  . 

°~U^ O)To(O)  [^(O)]2  *  1  ' 

H  I'.ever,  if  $2  is  to  make  a  proper  match  with  the  Y  =  0(1)  region,  we  must  have 

*1  =  AoM5(0)]P.  (8.42) 

It  turns  out  that  the  logarithms  are  crucial  for  the  generation  of  non-neutral  values  of 
Co.  The  solution  of  (8.39)  may  be  written 


=  A„[£7;(0)K-c1]!n 


c,  -  Ul(0)Y 


(8.43) 


(8.44) 


where  we  have  only  retained  terms  which  are  of  immediate  concern.  The  crucial  significance 
of  (8.44)  is  the  presence  of  the  -f -iir  jump  in  the  value  of  the  logarithm  (Mack,  1984,  for 
example),  as  the  point  F  =  Ci/Uq{0)  is  traversed;  it  is  this  that  generates  an  imaginary 
component  to  $2,  given  as  Y  00  by 


*5  ~  AoiTO0)F  -  co], 

(here  and  below  a  superscript  iV  denotes  an  imaginary  component). 

Returning  to  our  discussion  of  the  core  region,  we  see  we  must  have  that 

$i(0)  =  Ml)  =  o> 


(8.45) 


(8.46) 


and  consequently  t\  is  determined  by  the  expression  (8.14),  with  A  =  0.  On  account  of  the 
nature  of  Ri,&2,R3,v+,p+  we  expect  cx  to  be  a  real  quantity. 

We  are  now  in  a  position  to  specify  the  following  boundary  conditions 


u2(l)  =  0,  u2(0)  =  —  A0ci7r,  (8.47) 

and  thus  may  be  determined,  in  principle,  by  means  of  the  condition 

/W+  -f  Rip+)dy  =  v(0)Aocxt,  (8.48) 

Jo 

where  superscript  ‘r’  denotes  a  real  part. 

Unfortunately,  we  see  from  (8.41),  that  A0  is  dependent  upon  Tq(0)  and  Uq( 0),  and  since 
both  these  quantities  are  zero  in  the  case  of  an  adiabatic  lower  wall,  then  in  this  case  we 
must  consider  the  expressions  (8.31)  up  to  0(a3)  in  order  to  determine  a  value  for  the 
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imaginary  component  of  the  wavespeed  as  a  — ►  0;  however,  since  the  general  technique  is 
well  established  above,  we  do  not  carry  this  out  in  this  paper. 

The  key  result,  therefore,  is  then  that  Jm{co}  is  0(a2)  in  general  (but  0(&3)  in  the 
insulated  wail  case),  and  as  a  consequence  we  do  not  expect  to  obtain  non- neutral  values  of 
ci,  from  (8.14)  as  ai  — »  oo  (this  would  only  be  expected  in  a  study  of  the  0(e2)  terms  of 
the  c  expansion,  in  general,  and  the  0(e3)  terms  in  the  case  of  an  insulated  lower  wall).  The 
above  also  clearly  illustrates  how  cq  is  only  complex  for  values  of  a  >  «o>  since  for  a  <  a0 
the  ‘jump’  in  the  value  of  the  logarithm  is  not  present. 

9.  A  Change  of  Boundary  Conditions 

It  is  curious  that  the  (confined)  flow  profile  under  consideration  in  this  paper  has  such 
different  stability  characteristics  from  those  of  the  (external)  boundary  layer.  This  leads  us 
to  question  the  nature  of  these  fundamental  differences. 

One  obvious  candidate  for  investigation  is  the  effect  of  the  boundedness  of  the  domain 
(i.e.,  the  impermeability  condition  imposed  on  two  walls).  In  an  attempt  to  assess  this  effect, 
we  considered  the  problem  where  the  impermeability  constraint  on  the  upper  wall  (y  =  1),  is 
replaced  by  one  of  radiation,  while  retaining  the  same  basic  profile.  Specifically,  we  replace 
(5.3)  by 

«o(0)  =  0,  (9.1) 

(as  before),  but  with 

+  on  y  —  l,  (9.2) 

where 

v  =  ±a[l-Ml(l-coY}*.  (9.3) 

Here  the  sign  of  v  xa  chosen  to  ensure  that  Rl{ v}  >  0.  Indeed  this  amounts  to  considering 
the  piecewise  continuous  flow  which  consists  of  our  basic  Couette  flow  for  0  <  y  <  1,  and 
uniform  flow  ( Uo{y )  =  T0(y)  =  1)  over  1  <  y  <  oo.  This  may  be  shown  using  the  arguments 
used  by  Drazin  and  Reid  (1981)  who  considered  incompressible  piecewise  linear  velocity 
profiles  and  imposed  continuity  of  pressure  at  y  =  1. 

The  system  defined  by  (5.1)  -  (5.2),  (9.1)  -  (9.2)  is  solved  numerically  by  the  techniques 
used  in  obtaining  the  results  of  Section  6  (i.e.,  a  Runge-Kutta  shooting  scheme  in  conjunction 
with  a  Newton  iteration  scheme  to  iterate  on  Co). 

Results  for  one  example  were  obtained,  namely  for  M =  2,  with  adiabatic  conditions  on 
the  temperature  at  y  =  0.  Results  for  Rl{co}  are  shown  in  Figure  23.  Two  modes  were  found, 
the  first  (mode  I)  originates  at  a  =  0,  with  a  wavespeed  Co  =  1  —  1/Moo  —  0.5,  and  as  such 
is  typical  of  so-called  “first  modes”  observed  in  supersonic  boundary  **  er  stability  studies 
(e.g.,  Mack,  1987).  Mode  II,  on  the  other  hand,  originates  with  a  wavespeed  co  «  1.48, 
and  consequently  is  more  typical  of  the  bounded  flow  modes  observed  in  previous  sections 
of  this  paper.  Figure  24  shows  the  variation  of  7ra{co}  with  a.  Mode  I  is  unstable  over  the 
entire  range  of  a  shown,  while  mode  II  is  neutrally  stable  up  to  a  «  3.5,  at  which  point 
Rl{co}  drops  below  unity,  and  7m{co}  becomes  negative,  indicating  a  stable  mode.  We 
were  unable  to  find  other  modes  over  the  range  of  a  considered,  particularly  modes  with 
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RI{cq}  <  0.  Thus  it  appears  that  the  imposition  of  radiation  type  boundary  conditions 
results  in  a  hybrid  situation,  with  stability  characteristics  similar  to  those  found  in  both 
bounded  and  unbounded  flows. 

10.  Conclusions 

In  this  paper,  we  considered  the  linear  stability  of  compressible  plane  Couette  flow.  Our 
numerical  results  of  the  full  governing  equations  are  clearly  in  agreement  with  the  predictions 
of  Morawetz  (1952,  1954),  and  also  there  is  a  qualitative  agreement  with  the  work  of  Glatzel 
(1989),  although  the  latter  author  considered  a  model  basic  flow. 

The  main  thrust  of  this  paper  has  been  directed  towards  modes  which  have  been  pre¬ 
dominantly  inviscid  in  nature.  It  appears  that  these  are  the  most  important,  since  they  may 
become  unstable,  while  our  studies  indicate  that  viscosity  plays  a  generally  stabilizing  role, 
throughout. 

The  expansions  of  the  assumed  form  (4.11)  are  typical  of  solutions  which  are  predom¬ 
inantly  inviscid,  but  have  viscous  corrections  (see  also  Morawetz  1952).  In  Section  8  we 
studied  the  important  regimes  where  Co  — »  0  (or  cq  — ►  1)  in  which  the  first  two  terms  in 
the  series  become  comparable,  a  critical  layer  forms  clwse  to  one  of  the  walls,  and  the  vis¬ 
cous  layer  thickens  to  0(Re~3 ).  These  are  the  “exceptional”  cases  referred  to  by  Morawetz 
(1952). 

In  Section  9,  it  was  shown  how,  to  a  large  extent,  it  is  the  boundary  conditions  imposed 
on  disturbances,  that  determines  the  nature  of  the  stability  of  the  flow.  Indeed,  earlier  in  the 
paper,  in  Section  4  it  was  shown  in  the  study  of  the  G.I.P.  condition  that  it  is  important  to 
impose  correct  (i.e.,  radiation  type)  boundary  conditions  in  the  case  of  unbounded  flows  to 
avoid  the  (erroneous)  requirement  that  a  G.I.P.  is  necessary  for  the  existence  of  supersonic 
neutral  modes. 
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Figure  2:  Phase  velocity  spectrum  at  =  2,  a  =  0.1,  Re  =  7  x  105.  (a)  75  pts,  (b)  100 
pts,  (c)  125  pts,  (d)  150  pts. 


30 


(a)  75  pts,  (b)  100 


100  pis 
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150  pts 


Figure  5:  Phase  velocity  spectrum  at  =  5,  a  =  0.1,  Re  =  2  x  106.  (a)  100  pts,  (b)  125 
pts,  (c)  150  pts. 
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Figure  6:  Phase  velocity  spectrum  at  =  5,  a  =  0.1,  Re  =  5  x  106.  (a)  100  pts,  (b)  125 
pts,  (c)  150  pts. 
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Figure  9:  CiVRe  as  a  function  of  a  for  Mm  =  2  (mode  II),  and  the  two  Reynolds  numbers 
Re  =  7  x  105  and  Re=  =  1.4  x  106.  The  third  curve  is  7m{ci}  obtained  from  asymptotic 
theory. 


Figure  10:  Cr  as  a  function  of  a  for  M ^  =  2  (mode  II),  and  the  two  Reynolds  numbers  Re 
=  7  x  105  and  Re=  =  1.4  x  106. 
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/m{co}  x  105 


Figure  13:  /m{co}  as  a  function  of  a  for  Moo  =  2,  adiabatic  lower  wall,  mode  II. 
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Figure  14:  RI{cq}  as  a  function  of  a  for  =  5,  adiabatic  lo”  nr  wall,  mode  II. 


41 


Rt{  co} 


Figure  15:  i2/{co}  as  a  function  of  a  for  =  5,  adiabatic  lower  wall,  modes  IA,  IB,  Ic. 
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Figure  21:  Jm{ci}  as  a  function  of  a  for  —  2,  adiabatic  lower  wall,  modes  I  and  II. 


a 

Figure  22:  Im{ci}  as  a  function  of  a  for  =  5,  adiabatic  lower  wall,  modes  I  and  II. 
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